Theoretical examination of stress fields in Pb(Zr .5Ti .5)C)3 
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In this paper, we develop a rigorous formulation of the local stress field. This approach can be used 
in conjunction with any first-principles method to study stress fields in complex bonded systems. In 
particular we investigate the induced stress fields resulting from the homogeneous deformations of 
tetragonal PbTiC>3 and rhombohedral PbZrOs. As an extension of these findings we also compute 
the induced stress fields resulting from homogeneous deformation of the (100) and (111) orderings 
of Pb(Zro.5Tio.5)03. The stress-field responses in these four materials are compared and their 
piezoelectric responses are discussed. 
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I. INTRODUCTION 



When mechanical stress is applied to a piezoelec- 
tric material, microscopic atomic rearrangements occur 
which give rise to a change in the macroscopic electric 
polarization of the material. Conversely, application of 
a voltage across a piezoelectric material produces an in- 
ternal strain within it. In either case, it is the material's 
internal structural behavior (and therefore spontaneous 
polarization) under applied or induced stress that under- 
lies the piezoelectric response. 

The ability of piezoelectric materials to interconvert 
electrical and mechanical energy lies at the foundation 
of many electro-optic and electro-acoustic devices. The 
use of piezoelectric materials in these devices stems from 
the need to monitor the magnitude of induced or applied 
electrical response through the device. One of the clear- 
est examples of such an application is seen in the constant 
tunneling voltage mode of the scanning-tunneling micro- 
scope (STM)Jd The STM scanning needle is attached to 
a cylinder of piezoelectric ceramic. A voltage is applied 
across the gap between the needle tip and the scanned 
surface. If there are any structural variations in the 
scanned surface, this will cause a change in the voltage. 
Using a feedback loop, the piezoelectric ceramic is either 
compressed or expanded by altering the current passing 
through the ceramic so as to adjust the height of the 
scanning tip relative to the surface until the originally 
applied voltage is restored. Monitoring the height of the 
tip provides information on the structure of the scanned 
surface. Piezoelectric materials also play a vital role in 
electro-acoustic transducers. In these devices the piezo- 
electric material act as an interpreter for the incoming 
(or outgoing) sound wave and the outgoing (or incom- 
ing) electric signal. These types of devices have uses in 
underwater and medical ultrasonic imaging. 

In this paper, we present a first-principles investiga- 
tion of the distorted perovskite materials PbTi03 and 
PbZr03 at zero temperature and study the spatial varia- 



PbTiOs, PbZrOa, PZT, (100)-Pb(Zr .5Tio.5)0 3 , (111) 



tion of their stress-field responses to an externally ap- 
plied uniform strain. In addition, stress-field stud- 
ies of the (100) and (111) orderings of the solid solu- 
tion Pb(Zr .5Tio.5)0 3 (PZT) are also reported. Our 
selection of PZT is motivated by the large quantity 
of experimental studies describing the existence of a 
strong piezoelectric response in various composite PZT 
ceramics.! 2 ] In Section II we briefly describe the formal- 
ism for construction of the local stress fields. In Sec- 
tion III we present results for the computation of lo- 
cal stress fields induced by a uniform uniaxial deforma- 
tion. A discussion of the local stress fields produced 
in the simple perovskite crystals and the more com- 
plex PZT superlattices is provided in Section IV. We 
review the main conclusions of this study in Section V. 

II. STRESS-FIELD FORMALISM 



It is central to the study of any piezoelectric crys- 
tal to understand the material's internal response to an 
externally applied strain. Experimentally, a particular 
crystal's structural response to an applied strain can 
be measured using various diffraction techniques! 3 ] or di- 
rect measurement of the changes in the dimension of the. 
sample by electrical capacitance or optical interference.!!! 
In order to ascertain the effect strain has upon piezo- 
electric response, these techniques are paired with high- 
field measurement of strain hysteresis and polarization 
change. More recently, field-induced strain has bee 
measured using a displacement magnification technique. 
However because the formation and testing of many 
of these strained materials is complicated and possibly 
destructive,! it is advantageous to have a concise the- 
oretical method from which information concerning the 
microscopic response of a system to an external strain 
can be extracted. 

We consider the response of a system to a homogeneous 
long wavelength deformation (scaling transformations de- 
scribing pure dilation, strain or shear). For an interacting 
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system of atoms, the introduction of any of these uniform 
deformations can induce a force distribution on all the 
structural degrees of freedom of the unit cell. Within a 
harmonic theory, the induced atomic force distribution, 
F, contains all the relevant information about the redis- 
tribution of the external stress within the cell. Therefore, 
the starting point for the computation of the local stress 
fields is the calculation of the induced force distribution. 
This can be accomplished using local density functional 
theory. These theoretical methods have proven to be very 
successful for studying structural phenomena in a broad 
class of condensed phases. Below we outline our method 
to compute the stress field given the local force density. 
Once the stress field has been constructed, correlations 
between the elastic response and the structural features 
of the system can be made. A more detailed explanation 
of the method will be presented elsewhere! 

We begin with the principle of virtual work in the pres- 
ence of an induced force distribution F m . Any set of 
displacements of the nuclear coordinates u for a particu- 
lar interacting system produces a variation of the energy 
according to 
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where m represents the m-th ion of the interacting sys- 
tem. It is useful to convert the displacements and forces 
to continuous fields: 
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In the last line we have used lattice translational symme- 
try to transform the force and displacement distributions 
to a reciprocal space representation. The principle of vir- 
tual work can be recast in terms of the internal strains 
in the structure 



(3) 
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the symmetry of e n . Transforming this relationship into 
Fourier space gives the particularly convenient result 



(-5) 



g ■ r„ ■ u 



(5) 



Combining equation (5) for all 6 components of the 
strain tensor gives 



V-u (G 



or 
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Inserting this relationship into the last line of equation 
(2) and equating lines (2) and (3) yields 



(8) 
(9) 
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Equating the arguments of the summations in equa- 
tions (8) and (9) yields all 6 components of the local 
stress field in reciprocal space: 



a (-G)=f(-G\-V- 1 



By transforming according to 
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we can construct the stress distribution in direct space, 
and this result can be used to generate a map of the 
spatial distribution of the n-th stress field of the system. 

Vanderbilto has correctly pointed out that this for- 
malism only enables computation of the stress-field com- 
ponents which possess the periodicity of the unit cell. 
In addition the G=0 stress tensor C|»p,be computed by 
the approach of Nielsen and Martin.EnilJ Development of 
a method for the computation of the stress-field com- 
ponents which are uniform in one or two dimensions 
and of finite wavelength in the others is in progress. 



where u is the 6-component stress field tensor and e is 
the 6-component strain tensor. (Since only the contrac- 
tion of two rank-2 tensors is required in equation (3), the 
tensors are represented as length-6 vectors for simplicity.) 

The direct space components of the strain tensor can 
be directly related to the displacement field, 



e« (f ) = V • r„ • u(r) 



(4) 



where e„ is the n-th component of the strain tensor 
(n = 1,...,6) and T n is the 3x3 matrix determining 



III. RESULTS 

The first-principles calculations presented in iiiis paper 
are performed within density functional theoryjiil and the 
local density approximation (LDA) is used to describe the 
electron-electron interactions. For the solid-state calcu- 
lations, the single electron wave functions are expanded 
in a plane-wave basis using a cutoff energy of 50 Ry. 

To describe the electron-nuclear interaction, opti- 
mizftd, pseudopotentialsa in fully separable nonlocal 
forrr£3 are used.pWVn additional feature of our nonlocal 
pseudopotentialstU is their improved transferability over 
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a wide range of electronic configurations. We have been 
able to exploit the flexibility contained in the separation 
of the local and nondocal parts of the pseudopotential. 
By designing the form of the local potential so that the 
pseudo-eigenvalues and all-electron eigenvalues agree at 
an additional charge state, it is possible to improve the 
transferability of the potential across the charge states ly- 
ing between the original reference state and this second 
charge state. 

Due to the need for high accuracy when examining 
ferroelectric phenomena, semi-core shells are included in 
the generation of the pseudopotentials. We include as 
valence states the 3s and 3p for Ti and the 4s and 4p for 
Zr. The 5d shell is included for Pb. Furthermore, scalar 
relativistic effects ap, included in the generation of the 
Pb pseudopotential.Ej For each metal, a pseudopotential 
is constructed using a designed local potential with the 
addition of a square well within the core region. By do- 
ing so, we are able to achieve excellent transferability of 
the pseudopotential over a variety of charge and excited 
states. For each of the metals, excellent agreement of the 
pseudo-eigenvalues and total energy differences with the 
all-electron results is achieved for charge states of +4 to 
neutral. The oxygen pseudopotential is constructed using 
the s angular momentum channel as the local potential. 

Brillouin zone integrations for PbTiC>3 and PbZrC>3 
were done r-u^ing a 4 x 4 x 4 Monkhorst-Pack k- 
point mesh.E3 For tetragonal PbTiC>3 and rhombohc- 
dral PbZrC>3 this yields 6 and 10 fc-points, respectively 
in the irreducible wedge of the Brillouin zone. For the 
PZT crystals, the same fc-point mesh was folded into the 
smaller Brillouin zone, resulting in 3 irreducible fc-points 
for (100)-PZT, which was taken to have tetragonal sym- 
metry, and 5 irreducible /c-points for (lll)-PZT, which 
was taken to have rhombohedral symmetry It should 
be noted that to compute the nonuniform force distribu- 
tion resulting from the application of a uniform external 
stress, additional /c-points were needed due to the broken 
symmetries in the distorted structure. 

The calculations involving the rhombohedral PbZrC>3 
and (lll)-PZT deserve closer consideration. Crystallo- 
graphically, the zero-temperature form ofJ^bZrOs is or- 
thorhombic with 40 atoms per unit cellJlZTO In order 
to simplify our comparisons of the local stress fields, the 
structure of PbZrC>3 was fully relaxed within the rhombo- 
hedral symmetry of zero-temperature ferroelectric phase 
of PbZrC>3 with small Ti doping.0 Rhombohedral strains 
of the unit cell shape away from cubic were neglected, 
since they have been shown to be quite small (~ 0.1°).E2l 
In the case of (lll)-PZT, rhombohedral strains were also 
neglected to simplify analysis of the induced local stress 
fields. 

As part of the local stress-field calculations, complete 
structural relaxations of both internal coordinates and 
crystal lattice constants were completed for the PbTiC>3 
and PbZrC>3 distorted perovskite structures. We report 
our atomic positions and lattice constants for tetragonal 
PbTiC>3 and rhombohedral PbZrC>3 in Table I and com- 



TABLE I. Computed and experimental equilibrium lattice 
constants and atomic positions for tetragonal PbTiOs and 
rhombohedral PbZrOa. 

Present Theory Experiment 



PbTiO :i 



a(A) 3.870 3.862 a 3.905 b 

c/a 1.063 1.054 1.063 

as(Ti) 0.531 0.537 0.540 

x(Oi,0 2 ) 0.604 0.611 0.612 

x(0 3 ) 0.098 0.100 0.112 

PbZrOa 

o(A) 4.143 4.12 c 

x(Zr) 0.540 0.545 

x(Oi) 0.057 0.061 

y,z (Oi) 0.583 0.590 



"Reference [20]. 
b Reference[21]. 
c Reference [22]. 

pare these results with previous theoretical and experi- 
mental values where possible. The absence of experimen- 
tal values for the rhombohedral PbZrC>3 is explained by 
our simplification of the PbZrC>3 structure as described 
above. In the case of PblCiOs, our agreement to Pf^vi- 
ously reported theoreticalEil and experimental resultsEj is 
quite good and is within the error expected from density 
functional solid-state calculations done within LDA. For 
the rhombohedral PbZr03, our results agree quite well 
with the theoretical results of Singh.E3 (In Singh's work, 
the lattice constant of rhombohedral PbZrC>3 was not re- 
laxed. The theoretically determined lattice constant of 
the high temperature cubic perovskite was used instead.) 

Atomic and lattice relaxations were also performed for 
the PZT superlattices. The experimental lattice con- 
stants as well as theoretical lattice constants and relaxed 
atomic positions are contained in Table II. Experimental 
values for the (100)-PZT ceramic are taken for the tetrag- 
onal 50-50 batch composition PZT ceramic according to 
Jaffe et al. (ceramic 3 using the notation from reference 
[2]). The experimental values for the (lll)-PZT ceramic 
are taken for a rhombohedral PZT ceramic close to the 
50-50 batch composition (ceramic 5 using the notation 
from reference [2]). 

We studied the induced local stress fields by calcu- 
lating the internal force distribution induced by a uni- 
form external strain. The force distribution is obtained 
from first-principles density functional the ory w ithin the 
LDA and the Hellmann-Feynman theorem E-ffia Since our 
stress-field formalism relies on the fact that any defor- 
mation must not take the system beyond lowest order in 
gradients of the total energy, attention must be paid to 
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TABLE II. Computed equilibrium lattice constants and 
atomic positions for (lOO)-PZT and (lll)-PZT. Experimental 
lattice constants are given for randomly ordered PZT ceramics 
close to the 50-50 batch composition. See text for description. 

Present Experiment 



(100)-Pb(Zr .5Tio.5)0 3 



o(A) 8.313 8.279 a 

c/a 0.480 0.487 

x(Pb 2 ) 0.468 

x(Ti) 0.211 

x(Zr) 0.714 

k(Oi) -0.049 

x(0 2 ,0 3 ) 0.185 

x(0 4 ) 0.424 

x(0 5 ,0 6 ) 0.660 

(lll)-Pb(Zro.5Tio.5)0 3 

a(A) 8.043 8.164 3 

x(Pb 2 ) 0.498 

x(Ti) 0.237 

x(Zr) 0.737 

x(Oi) -0.014 

z(Oi) 0.221 

x(0 4 ) 0.470 

z(0 4 ) 0.721 



a Reference[2]. 

the magnitude of the deformations. We have found that 
±0.2% deformations in lattice lengths and ±0.5° in lat- 
tice angles are within the harmonic limit of the potential 
energy. Extension of this work beyond harmonic order is 
a promising direction for the study of recently discovered 
piezoelectric single crystals which exhibit large reversible 
strains E3 

There are 6 homogeneous deformations which can be 
made to any crystal: dilation having the symmetry of 
x 2 + y 2 + z 2 , uniaxial strains — tetragonal with symmetry 
2x 2 — y 2 — z 2 and orthorhombic with y 2 — z 2 symme- 
try, and elementary shear operations (xy, xz, and yz). 
Based on crystal symmetry considerations, certain ho- 
mogeneous deformations are degenerate and can easily 
be constructed from the other deformations. For a crys- 
tal subjected to any of the 6 homogeneous deformations 
one finds a local internal stress field in all 6 stress com- 
ponents. It is important to note that for a given deforma- 
tion, the induced stress fields corresponding to the 5 other 
deformations must each integrate to zero over the entire 
unit cell. As an example, for a tetragonal deformation, 
there is, of course, a tetragonal local stress-field response. 
In addition, one induces local dilation, orthorhombic and 
shear responses. The latter 5 responses, due to symmetry 
constraints, must each integrate to zero. 



For brevity, we only report the local stress fields in- 
duced by an applied tetragonal uniaxial stress. In partic- 
ular we focus our discussion here on the internal dilation 
and tetragonal stress-field responses induced by a tetrag- 
onal compression. A more comprehensive presentation of 
all the induced local stress fields originating from each of 
the 6 deformations is in progress E3 

To simplify the visualization of the stress fields, we 
have chosen to show only the regions of highest induced 
local stress. For the dilation response to a uniaxial com- 
pressive deformation, light blue isosurfaces correspond to 
regions in the unit cell undergoing compression in all di- 
rections and pink isosurfaces correspond to expansion in 
all directions. For the tetragonal response to a tetragonal 
compressive deformation, light blue regions correspond to 
prolate response (expansive along the axial direction but 
compressive along the equatorial directions) and pink re- 
gions correspond to oblate response (compressive along 
the axial direction but expansive along the equatorial di- 
rections) . 

A. Tetragonal PbTi0 3 

Figures IA and IB show the dilation and tetragonal 
local stress fields, respectively, produced in response to a 
uniaxial tetragonal deformation of the zero-temperature 
equilibrium structure of tetragonal PbTiC>3. The (100) 
lattice direction contains the Pb atoms at the lower left 
forward corner and upper left forward corners of the 
unit cell. It is instructive to recognize the regions under 
the highest induced stress for both stress components lie 
along the (100) direction, as this coincides with the di- 
rection of ferroelectric distortion in the relaxed crystal 
and the direction of applied stress. 

Dilation Response: The most salient feature of this 
response is its location relative to the Pb and Ti atoms. 
We find the region of highest induced stress does not in- 
volve the Pb atoms but instead straddles the Ti atom 
and is oriented along the (100) direction. The spatial 
representation of this component of the local stress field 
indicates that under uniaxial tetragonal compression, the 
response of the crystal would be to compress the volume 
in the half of the TiC>6 octahedron lying above the Ti 
atom. A commensurate expansion of the volume in the 
half lying below the Ti atom is also seen in this response. 
The combination of these volume deformations would re- 
sult in the shifting of the Ti atom toward the upper half of 
the octahedron. In contrast, the oxygens lying collincar 
with the regions of highest induced stress would shift in 
the (100) direction, opposing the motion of the Ti atom. 
This oxygen motion is also indicative of an expansion of 
the lower half of the TiOg octahedron and a compression 
of the upper half. Finally, there is no significant stress- 
field response involving the oxygens lying equatorial to 
the Ti atom. 

Tetragonal Response: Once again we find the region 
of highest induced stress does not involve the Pb atoms 
but surrounds the Ti atom and lies along the (100) direc- 
tion. Although the interpretation of this stress compo- 
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nent is different from the dilation response, its net atomic 
response is similar. In the upper half of the TiC>6 octa- 
hedron lying along the (100) direction, there is an oblate 
response which would shift the Ti along the (100) direc- 
tion toward the upper half of the octahedron. In addition 
to the oblate response, there is a prolate response in the 
lower half of the Ti06 octahedron. The effect of these 
volume deformations would be to elongate the lower half 
of the oxygen octahedron while shortening the upper half. 
The resulting oxygen motion is consistent with the vol- 
ume deformations. 

B. Rhombohedral PbZr0 3 

Figures IIA and IIB show the dilation and tetrago- 
nal local stress fields produced in response to a uniax- 
ial tetragonal deformation of the zero-temperature equi- 
librium structure of rhombohedral PbZrC>3. The (111) 
atomic distortion direction contains the Pb atoms at the 
lower left forward and upper right rear corners of the unit 
cell. It is useful to recognize that some regions under the 
highest induced stress for both stress components lie not 
exclusively along the (111) distortion direction, but also 
along the (100) direction, the direction of uniaxial strain. 

Dilation Response: In this stress field, as in the dila- 
tion field for PbTiC>3, there is pairing of compressive and 
expansive regions surrounding the central metal atom. 
However, there is a clear difference between the dilation 
response for the 2 crystals. In PbZrC>3 response we find 
a compression of the volume above the Zr and an expan- 
sion of the volume below the Zr atom. These regions of 
the stress fields are some combination of the (111) fer- 
roelectric distortion direction and (100) uniaxial strain 
direction. The corresponding motion of the Zr is toward 
the upper half of the ZrOg octahedron, along a direction 
lying between the (111) and (100) directions. In addi- 
tion to the motion of the Zr atom, the O atoms lying 
axial to the Zr atom are moving in a direction opposing 
the Zr atom motion. This is clearly seen in the combi- 
nation of the compressive and expansive deformations of 
the volumes within the unit cell. 

Tetragonal Response: This response has some fea- 
tures in common with the dilation response of the 
PbZrC>3 crystal described above. Here again there is a 
pairing of regions oriented along a combination of the 
(111) and the (100) directions. These regions are also 
centered at the Zr atom. Specifically, we find a prolate 
region in the lower half of the ZrC>6 octahedron and an 
oblate region in the upper half of the octahedron. The 
resulting motion of the O atoms lying axial to the Zr 
atoms is similar to the the motion of the O atom in the 
dilation response of this crystal. 

C. (100)-Pb(Zro. 5 Tio. 5 )03 

Figures IIIA and IIIB show the dilation and tetrago- 
nal local stress fields produced in response to a uniax- 
ial tetragonal deformation of the zero-temperature equi- 
librium structure of (100)-Pb(Zr .5Tio.5)0 3 . The (100) 
lattice direction contains the Pb atoms at the lower 



left forward corner and upper left forward corners of 
the unit cell. (Our restriction to tetragonal symmetry 
forces the ferroelectric distortion to be along (100).) In 
both reported responses, the highest induced stress field 
lies along the (100) direction parallel to the ferroelectric 
distortion direction and the (100) direction of uniaxial 
strain. 

Dilation Response: For this super-cell, an alternat- 
ing pattern of compressive and expansive responses is 
found lying along the 4-fold rotation axis. In this re- 
sponse to the tetragonal uniaxial compression, the high- 
est induced stress fields are reminiscent of the stress field 
in the dilation responses of both distorted perovskite 
structures. However, while the response in the PbTiC>3 
sub-unit is remarkably similar to the dilation response of 
the tetragonal PbTi0 3 described above, the response in 
the PbZrC>3 half of the unit cell shows an important dif- 
ference from the dilation response of the rhombohedral 
PbZrC>3. In the simple distorted zirconate, Zr atomic mo- 
tions along a composite direction lying between the (100) 
and (111) directions became active under tetragonal de- 
formation. As a result of our imposition of tetragonal 
symmetry, we find response in the (100)-PZT superlat- 
tice lying exclusively along the (100) direction indicating 
atomic motions purely in the (100) direction. In each 
oxygen octahedron, there is a compression of the lower 
half along with an expansion of the upper half. This leads 
to Ti and Zr motions toward the lower halves of the octa- 
hedra. The axial O atoms show opposing motion to the 
Ti and Zr atoms. 

Tetragonal Response: In this response to the uni- 
axial deformation, the induced stress field shows an al- 
ternating pattern of prolate and oblate components lying 
along the 4-fold rotation axis. The resulting motion of 
the atoms is consistent with the dilation response. 

D. (lll)-Pb(Zr .5Tio.5)0 3 

Figures IVA and IVB show the dilation and tetrago- 
nal local stress fields produced in response to a uniax- 
ial tetragonal deformation of the zero-temperature equi- 
librium structure of (lll)-Pb(Zro. 5 Tio.5)0 3 . The (111) 
atomic distortion direction contains the Pb atoms at the 
lower left forward and upper right rear corners of the unit 
cell. The (100) lattice direction contains the Pb atom at 
the lower left forward and upper left forward Pb atoms 
in the unit cell. There are 2 crystallographically unique 
Pb atoms lying along the (111) direction in addition to 
the Ti and Zr atoms. Due to ferroelectric distortion, the 
Ti atom has shifted along (111) toward the Pb atom at 
the lower left forward corner and away from the Pb in 
the center of the figure. The Zr atom has moved toward 
the central Pb atom along the (111). 

Dilation Response: Each individual region of high- 
est induced stress field in the response is oriented perpen- 
dicular to the (111) atomic distortion direction. Unlike 
the responses for the previously mentioned crystals, we 
find the Ti and Zr atoms of (lll)-PZT have negligible in- 
volvement. Instead we report alternating expansion and 
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compression regions consisting of entire PD2O4 octahe- 
dra. These octahedra are comprised of 2 adjacent Pb 
atoms lying along the (100) lattice direction and 4 equa- 
torial O atoms. The overall pattern of the expansions 
and compressions shows the (111) stacking of the entire 
superlattice. Focusing on the Pb at the lower left for- 
ward corner, we find expansion of the Pb2C>4 octahedron 
above this atom and contraction of the Pb2C>4 octahe- 
dron below, resulting in largely downward motion of this 
Pb atom. Analogous reasoning shows that the dilation 
stress field causes upward motion of the central Pb atom. 
The expansion and contraction of the Pb2C>4 octahedra 
also gives rise to equatorial oxygen atomic motion in the 
(100) plane. The O atoms not involved in the Pb2C>4 
octahedra also exhibit atomic motion. Due to the ellip- 
soidal shape of the stress-field regions, the O atoms lying 
beneath the Ti atoms are moving primarily upward along 
the (100) while the O atoms lying beneath the Zr atoms 
are moving mostly downward. 

Tetragonal Response: In this response we report a 
feature of the stress field not found in any of the other re- 
sponses for the other materials. Namely, we find prolate 
regions in the titanate sub-units alternating with oblate 
regions in the zirconate sub-units. In particular, the re- 
sponsive regions encompass the entire MC>6 octahedra of 
their respective sub-units. In each Zr sub-unit, there is 
a prolate response while in the TiC>6 octahedron there is 
an oblate response. 

IV. DISCUSSIONS 

As a general guideline for analysis, it is important to 
understand the relationship between the ferroelectric dis- 
tortion for a particular material and the direction of uni- 
axial stress. It is the combination of these effects that 
will either enhance or diminish the effect of the ferroelec- 
tric distortion and therefore the extent of piezoelectric 
response. For all of the above reported responses the 
direction of uniaxial stress is (100). 

As previously stated, the ferroelectric distortion for 
PbTi0 3 lies along the (100). Although pure PbTi0 3 
is not piezoelectric, piezoelectric response has been re- 
ported for a PbTiC>3 ceramic doped with 2.5 mol% 
LaOo.5 and 1.0 mol% M11O2.U While the piezoelectric re- 
sponse of this material is smaller than that of the PZT 
solid solution ceramics, it is interesting to note that the 
calculated stress fields for the undoped material show 
the ability of the crystal to exhibit such a response upon 
application of a uniaxial strain. By examination of the 
results in Table I, it is clear that the relative positions 
of the Ti and O atoms underlie the ferroelectric response 
in this crystal. In particular, we find the Ti atom lies 
~0.28A below the equatorial oxygen plane. This means 
that in the zero-temperature tetragonal phase the vol- 
ume beneath the Ti atom is compressed while the vol- 
ume above the Ti atom is expanded with respect to the 
cubic prototype structure. We have found in both re- 
ported responses a tendency for the volume beneath the 



Ti atom to expand and the volume above the Ti atom 
to compress upon application of a uniaxial compression. 
These volume deformations are indicative of a structural 
phase transformation from tetragonal to a para-electric 
cubic phase. Since the cubic structure shows no ferro- 
electric response, the application of a strain that returns 
the tetragonal PbTiC>3 toward the cubic phase will re- 
duce the ferroelectric properties and therefore underlies 
the piezoelectric response of this material upon uniaxial 
compression. 

The structure of rhombohedral PbZrC>3 shows an 
atomic distortion along the (111) direction. Upon ap- 
plication of the tetragonal strain, we find the highest in- 
duced stress field lying along a combination of the strain 
direction and the atomic distortion direction. It is clear 
in these responses that, for this crystal, the directions 
of uniaxial strain and atomic distortion are competitive 
influences on the atomic displacements. Despite the com- 
plexity of the predicted atomic motions, the overall re- 
sponse displays the existence of a piezoelectric response 
in this material. The manifestation of the piezoelectric 
response may be a result of the relative closeness in en- 
ergy of the anti- ferroelectric ground state of PbZrOaES 

For the (100)-PZT ceramic, we report an atomic dis- 
tortion along the (100) direction. As stated in Table II, 
the Ti atom is positioned ~0.22A above the plane con- 
taining the equatorial O atoms. Also the Zr atom lies 
~0.45A above its equatorial O atom plane. This means 
that before application of the tetragonal strain, the up- 
per halves of the MOg octahedra are compressed and 
their lower halves expanded. Upon application of com- 
pressive uniaxial strain, we again find both expansion of 
the halves of the octahedra containing metal atoms and 
compression of the empty halves. This results in metallic 
and oxygen motion opposing the ferroelectric distortion. 
These atomic motions are very similar in nature to the 
reported responses in the tetragonal PbTiC>3. Once again 
we note that the application of the uniaxial compression 
is done in a manner that will lead to a structural trans- 
formation towards a ceramic comprised of cubic-like sub- 
units. A deformation that favors a para-electric phase 
will decrease the ferroelectric effect. Therefore we con- 
clude that application of a uniaxial tetragonal compres- 
sion causes a piezoelectric response in (100)-PZT. It is 
important to note that the structure of our (100)-PZT 
maintains the 4-fold rotation axis, a feature that would 
not allow any atomic motions along (111). We may find 
composite atomic motions in the zirconate sub-unit, sim- 
ilar to the motions predicted in the PbZrOs crystal, upon 
removal of this symmetry constraint. 

The atomic distortions for the rhombohedral (111)- 
PZT equilibrium structure lie along the (111) direction. 
In the undeformed superlattice, we find significant Ti 
and Zr displacements along the (111) relative to the Pb 
atoms. Upon uniaxial compression we do not find any 
motion of the Ti or Zr atoms. Instead we find complex 
Pb and O motions. We report alternating motions of Pb 
atoms primarily along (100) as well as O motions along all 
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three lattice directions. The origin of O motion in (111)- 
PZT deserves particular attention. The motion along 
(100) and (100) of the oxygen atoms that lie directly 
above and below the Zr and Ti atoms can be predicted 
from either the dilation or tetragonal responses. How- 
ever, the complex motion of the other oxygen atoms can 
only be ascertained from a simultaneous analysis of both 
responses. As an illustration, we examine the motion of 
the O atom located near (|,0,|). Because of its position 
relative to the expansion region depicted in the dilation 
response, it has a component of atomic motion along the 
(001) direction. However, according to the tetragonal 
response, the same atom also moves along the (010) di- 
rection due to the pairing of an oblate response centered 
in the TiOg octahedron and a prolate response centered 
in the adjacent ZrOg octahedron. This same type of anal- 
ysis can be done for all the oxygen atoms lying equato- 
rial to the Zr and Ti atoms and can be used to predict 
their atomic motions due to uniaxial compression. Based 
upon these atomic motions, it is possible to characterize 
the piezoelectric response in (lll)-PZT. The response 
can be summarized as complex motions of the Pb and O 
sub-lattices against the almost stationary Zr/Ti atomic 
positions, and not a simple move toward a higher symme- 
try structure, as seen in the previously described crystals. 

V. CONCLUSIONS 

In this density-functional study, we have examined the 
induced local stress fields resulting from an externally ap- 
plied homogeneous uniaxial deformation of two distorted 
perovskite structures, PbTiOs and PbZrOa. We have 
found in the case of tetragonal PbTiOs the existence of a 
piezoelectric response involving atomic motions of the Ti 
and axial O atoms upon application of a tetragonal com- 
pressive strain. We also report a piezoelectric response in 
the rhombohedral PbZrOa due to atomic motions along 
the (111) atomic distortion and (100) uniaxial strain di- 
rections. In addition we have investigated the stress fields 
of the (100)-PZT and (lll)-PZT superlattices. We find 
that the stress responses of the (100) superlattice dis- 
play similar features to those of the simpler distorted 
perovskites, most notably in the motion of the Ti and 
Zr atoms. However in the (lll)-PZT crystal, we have 
identified features not represented in any other crystal 
studied. We find complex Pb and O motions against a 
fixed Zr/Ti sub-lattice upon application of a tetragonal 
stress. Application of stress along a direction not paral- 
lel to the superlattice ordering and ferroelectric distortion 
directions gives rise to a complicated stress-field pattern. 

This study demonstrates the utility of stress fields 
and their applicability to ferroelectric phenomena. The 
construction of the local stress field is computationally 
tractable and provides an intuitive way to visualize and 
understand the response of a crystal to applied stress. 
Furthermore, we feel that stress-field analyses done in 
conjunction with spontaneous polarization studies will 
broaden the understanding of piezoelectric materials. 
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FIGURES 

FIGURE IA. Dilation stress-field response to a uni- 
axial tetragonal deformation of tetragonal PbTiC>3. 
(blue=lead, red=titanium, green=oxygen) . 

FIGURE IB. Tetragonal stress- field response to a uni- 
axial tetragonal deformation of tetragonal PbTi03. 
(blue=lead, red=titanium, green=oxygen) . 

FIGURE IIA. Dilation stress-field response to a uni- 
axial tetragonal deformation of rhombohedral PbZr03. 
(blue=lead, yellow=zirconium, green=oxygen) . 

FIGURE IIB. Tetragonal stress- field response to a uni- 
axial tetragonal deformation of rhombohedral PbZrC>3. 
(blue=lead, yellow=zirconium, green=oxygen) . 

FIGURE IIIA. Dilation stress-field response to a uni- 
axial tetragonal deformation of (100)-Pb(Zr . 5^0.5)0-3. 
(blue=lead, yellow=zirconium, red=titanium, green= 
oxygen). 

FIGURE IIIB. Tetragonal stress-field response to a uni- 
axial tetragonal deformation of (lOO)-Pb(Zro. 5^0.5)03. 
(blue=lead, yellow=zirconium, red=titanium, green= 
oxygen). 

FIGURE IVA. Dilation stress-field response to a uni- 
axial tetragonal deformation of (lll)-Pb(Zr . 5^0.5)0-3. 
(blue=lead, yellow^zirconium, red=titanium, green= 
oxygen). 

FIGURE IVB. Tetragonal stress-field response to a uni- 
axial tetragonal deformation of (lll)-Pb(Zro. 5^0.5)03. 
(blue=lead, yellow=zirconium, red=titanium, green= 
oxygen). 
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